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We explore the possibility of studying B nlv semileptonic decays at large recoil momentum. 
Our methods include the use of a random-wall source for the pion to reduce statistical errors, and 
different smearing functions are used for the B meson to improve the overlap with the ground 
state. We observe, in general, a factor of 3-4 improvement in the signal-to-noise ratio in correla- 
tion functions if random-wall propagators are used. 
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Figure 1: Kinematics of B ^ nlv semileptonic decay, ps is the momentum of the B meson and pn is the 
momentum of the pion, q is the momentum transfer. 

1. Introduction 

Precise determination of the form factors /o(^^), f+{q^) in B nlv semileptonic decays is 
crucial to the determination of the CKM-matrix element \Vub\- Fig. 1 shows the kinematics of 
the process. Lattice QCD provides a first principles nonperturbative approach to calculate the 
from factors in semileptonic decays. Standard simulation methods, however, are problematic in 
the low region when the pion has large recoil momentum. Simulation results are limited to 
> 15GeV^ ^ while experimental data spans the entire q^ range. Large recoil momenta are 
difficult for lattice calculations because statistical errors, which are set by E{p) — E{p = 0), become 
worse when the hadrons have large momenta. In addition, discretization errors, which are set by 
a^p^, increase as the pion momentum increases. 

In order to utilize all the experimental data and thereby reduce the experimental error on \Vub\, 
it is important to develop new simulation techniques to cover the low q^ region. One method is to 
reduce the pion momentum by using a lattice frame in which the B meson is moving in the opposite 
direction to the pion; to describe a b quark with large velocity on the lattice the "moving-NRQCD" 
formalism is used 0, ||, |6|]. In this work we concentrate on reducing the statistical errors with the 
use of a random-wall source for the Ught quark propagator. In particular we explore the possibility 
of simulating at cipji — ^(3,0,0), which corresponds to q^ ^ lOGeV^ in the B meson rest frame. 
We test this on MILC coarse lattices where discretization errors are about 8% at this momentum; 
errors will be much smaller on the fine (3%) and super-fine (1.5%) lattices. We also use the Highly 
Improved Staggered Quark (HISQ) action [^, ^ for the valence light quarks. 

2. Pion 2-Point Function 

A zero-momentum random-wall source is generated by setting the color vector on each site of 
a time slice' to a three-component random complex unit vector fi{x) [^]. This random source is 
used in matrix inversion to obtain the staggered-quark propagator g{y,x) 



(2.1) 



X X 



In practice we also put the source on a random time slice. 
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where My^x is the kernel of the staggered-quark action. It is more convenient to work in the naive- 
quark basis. The 4-component naive-quark propagator S{y,x) is given by ^ 

S{y,x) = (v/'(y)TF'W> =g{y,x)a{y)aUx), (2.2) 

where 

^{x) = ll{Y,yr (2.3) 

^1=0 

The pion correlator in the naive-quark basis becomes 

</5(y)/5W> = ((v7'(y)75V(3'))(yW75'/^'W)> 
= Tr[Y5w'{y)W'{x)Y5¥'ix)W'{y)] 
= Tr[Y5S{y,x)Y5Y5SHy,x)Y5] 
= Tr [£l{y)a^ {x)a{x)a^ {y)] tr [ \g{y,x)\^] 
= 4tr[\g{y,x)\^] (2.4) 

in terms of the staggered-quark propagator. Here Tr[. . .] is a trace over spinor indices while tr [. . .] is 
taken over the color indices. Therefore to construct the pion correlation function at zero momentum 
we take the magnitude square of g{y) at the sink and sum over spatial sites, and divide by the 
number of sites N 

^Ls*iy)s{y) = ^ E g*{y,x')g{y,x)^^*{x')^^{x). (2.5) 

y y,x,x' 

Since the averaged correlator has contributions only from where the quark and antiquark start at 
the same spatial site, i.e., x = x', the random-wall source simulates many-point source, and thereby 
increases the statistics. 

For correlator at finite momentum k, an additional phase is added to the source 

t{y) ^ l^g{y,x)e^''^'ri{x) = EM-,ie±'^T7(x). (2.6) 

X X 

In this case we multiple g^{y) and {g^(y))* with an explicit insertion of e'*'^ at the sink and sum 
over spatial sites. 

Fig. 2 compares the pion 2-point functions obtained with a local and a random source at 
apjt = (0,0,0) and apji = ^(3,0,0) on MILC coarse lattices. The advantage of using a random- 
wall source is clearly demonstrated, with statistical errors about 5 times smaller compared to the 
local case. To fit the 2-point functions we employ a Bayesian technique [ p^ and do a multiple- 
exponential fit 

d'\t)=W-lfate-^'^'. (2.7) 

k=Q 

We use 5 exponentials Njc = 5 (i.e., 3 normal states and 2 oscillating states). We observe a factor 
of 2-3 improvement in the ground state energy and the amplitude for the random-wall source. 
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Figure 2: Pion 2-point functions obtained with a local source and a random source on MILC coarse lattices 
[am^ea = -01 /.05, am^ai = .05465 (~ rris) for apjt = (0,0,0); am^ea — -02/. 05, am^ai = .02675 {'^ m^jT) for 
apn — Pi'- results for ground state energy £0 and amplitude ao are also shown. 



3. Heavy-Light 2-Point Function 

The heavy-light coirelator is 

= Tr[Y5Y'{x)¥'{y)75¥^{y)¥^{^)] 
= Tr[Y5Y5S'^{y,x)Y5Y5G{y,x)] 

= Tr[ {a\y)g*{y,x)) G{y,x)a{x)] . (3.1) 
Here G{y,x) is the heavy-quark propagator. We use NRQCD for the b quark 

1) . (1 - ^) (, - ^)"./(.) (■ - g)7, - f ) G(.,,). (3.2, 



and H = Ho + dn is the 0{\/M) improved lattice NRQCD Hamiltonian [[U], yi- In order to 
combine with the light-quark propagator, the heavy-quark propagator must be initialized with the 
same random noise ?] {x) 

G{x,t = 0) = Y,^{x'M\x'-x\)ri{x'), (3.3) 

y 

where —x\) is the smearing function centered atx'. Since the smearing function depends on 
the magnitude of x' — x only, G{x,t = 0) can be computed efficiently by applying translation to 
<I>(|x|), the smearing function with center located at the origin, when doing the summation. At the 
sink we multiply the heavy-quark propagator with D.\y)g*{y) x^{\y\) and sum over spatial sites 
and divide by A'^. Again only the terms with matching random numbers survive after the ensemble 
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Figure 3: Heavy-light 2pt-functions on MILC coarse lattices [amsea = .01/. 05, am^d = .05465 (~ mj), 
arrih = 2.8]. Gaussian smearing with radius ar = 4 is used and results are shown for local-source, local-sink 
[LL] and smeared-source, smeared-sink [SS]. A matrix fit of LL, LS, SL, SS gives [local] Eq — 0.5532(15), 
ao = 0.1402(10), and [RW] Eq = 0.5559(7), oq = 0.1424(5). 



average, and the correlation function is the average of all contributions from where the heavy quark 
and light quark start at the same site. 

Fig. 3 shows the results obtained with a local source and a random source, with and without 
smearing. The random-wall results are again more accurate, although the improvement is not as 
significant as for the light-light correlators. Gaussian smearing is used and results are shown for 
local-source, local-sink [LL] (i.e., no smearing) and smeared-source, smeared-sink [SS]. The plots 
show that the advantage of using random-wall propagators decreases if smearing is used. It is 
therefore crucial to fit the local and smeared results together. We fit LL, LS, SL, SS simultaneously 
in a matrix fit using 5 exponentials (3 normal states and 2 oscillating states) and find that a smearing 
function of radius ar = 4 gives the least relative errors for the ground state energy and amplitude. 



4. 3-Point Function 



We are interested in the 3-point function 

= Tr [75 V {y)W' (z) {z)W^ {x) 75 V i^W )] 
= Tr[75755t(z,3')757MG(z,x)755(x,3;)] 

= Tx[{a\z)g*{z,y))y5yMz,x)ys{^{x)g{x,y))]. (4.1) 

with xq < zq <yo- To compute the current we initialize G{z,x) with 75 x Q.{x)g{x) and then 
propagate the heavy quark backward^ in time to the current insertion point z, where it turns into a 



Backward means to the opposite direction of the light-quark propagator. 
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Figure 4: FIG. 5. The 3-point functions (corresponding to the temporal vector currents) obtained with 
local and random- wall propagators on MILC coarse lattices [the same ensembles as in Fig. 2]. The B 
meson starts at f = and the pion is at f = 16. The random-wall results are 3-4 more accurate compared 
to the results obtained using local propagators. In particular we obtain, for no smearing, v'oo = 0.0605(39) 
[local], voo = 0.0597(7) [RW] for ap,, = (0,0,0), and voo = 0.029(49) [local], vqo = 0.043(33) [RW] for 
= f (3,0,0). 



light quark; more precisely, we take the trace of the product of G{z,x) and n^(z)g*(z) x 757^ at the 
insertion point z- This setup is not merely convenient and efficient, but the same computer code can 
be used for local and random-wall propagators^ . We only have to divide the random-wall results 



by an extra factor of N since there are N times more contributions like the one in Eq. ( |4. 1| ), each 
from a different source point x, to the averaged correlator. 

In Fig. 4 we plot the temporal vector currents calculated with local and random-wall prop- 
agators for pion momenta ap^ = (0,0,0) and ap^ = ^(3,0,0). Results again clearly show that 
statistical noises can be suppressed substantially by the use of random sources. The 3-point func- 
tion has the functional form 



C(3)(0= £ £ (-l)^^(-l)'(^-)a,w 
k=o 1=0 



-E'it -E>,{T-t) 



(4.2) 



where T =yQ—xo. To extract the matrix element vqo we fit the 3-point function, the 2-point function 
(Fig. 2) and the heavy-light 2-point function (Fig. 3) simultaneously using Nji = Nb = 5. Fit results 
are given along with the graphs. We find that random-wall results are about 5 times more accurate 
compared to local case at apji = (0,0,0), and about a factor of 2 better at apji = ^(3,0,0). 



^In the local case we read in local propagators g(x) = £,,/ g{x,y')5yY 
case we read in the random- source propagators Eq. (p. lb. 



--g{x,y) and; 



- si^^y)' in the random- wall 



6 



B Semileptonic Decays at High Recoil Momentum 



K. Y. Wong 



5. Conclusion 

The major problem in studying semileptonic decays on the lattice is the exponential growth 
of statistical errors in correlation functions as the pion momentum increases. In this work we ex- 
plored the possibility of reducing statistical noises using random-wall light-quark propagators. We 
demonstrated that correlation functions obtained with random sources have much better signal-to- 
noise ratios, with statistical errors about 3-4 times smaller than those obtained with local propa- 
gators. One should therefore improve the current lattice simulations at small pion momenta using 
random-wall propagators. Encouraging results were also obtained at pion momentum as large as 
apji = ^(3,0,0) (corresponding to ~ lOGeV^), although statistical errors are still too large for 
the results to be useful. Work is in progress, e.g. fitting correlation functions with different T 
simultaneously, to further improve the calculations. 
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